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"q' Abstract 

This paper applies meshless method of lines, which uses radial basis functions (RBFs) as a spatial collocation scheme to solve the 
Coupled Drinfeld's-Sokolov-Wilson System. Runge-Kutta method is used for time integration of the system of ODEs obtained as 

_ a result of spatial discretization in contrast to usual RBFs or finite difference methods. Accuracy (L2 and L M ) is compared with the 

^"'existing results from other methods available in the literature. 
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1. Introduction 

Nonlinear partial differential equations model many impor- 
tant physical, chemical and biological phenomena, and their use 
recently has spread into economics, finance, image processing, 
medicine and other fields. In order to investigate the prediction 
of these models, it is often necessary to approximate their so- 
lution numerically. Numerous numerical methods are available 
in the literature; the mesh based methods including, finite dif- 
ference (FDM), finite element(FEM) and finite volume method 
(FVM), are most widely used, however, the construction of an 
appropriate mesh in arbitrary geometry is hard and expensive. 

Consequently, to avoid the mesh generation, in recent years, 
meshless techniques have attracted the attention of researchers 
as alternatives to traditional finite element, finite volume and fi- 
nite difference methods. In a meshless (meshfree) method a set 
of scattered nodes, with no connectivity information required 
among the set of points, is used instead of meshing the domain 
of the problem. Examples of some meshless schemes are the 
elementfree Galerkin method, the reproducing kernel particle, 
the local point interpolation, etc (e.g. see [? ] and references 
therein). 

Over the last two decades, the radial basis functions method 
is known as a powerful tool for scattered data interpolation 
problems. The use of radial basis functions as a meshless 
procedure for numerical solution of PDEs is based on the col- 
location scheme. Due to the collocation technique, this method 
does not need to evaluate any integral. The main advantage of 
numerical procedures that use radial basis functions over tra- 
ditional techniques is the meshless property of these methods. 
Radial basis functions are actively used for solving PDEs ([? 
0], and reference therein). In the above cited work, RBFs are 
used to replace the function and its spatial derivatives, while 
finite difference scheme is used to march in time. This method 
was first introduced by Kansa [? ? ] in 1990 for the numerical 



solutions of the PDEs. Kansa used the Multiquadric (MQ) RBF 
to solve the elliptic and parabolic PDEs. Recently, Flyer and 
Wright [? ] indicated RBFs allowed for a much lower spatial 
resolution, while being able to take unusually large time-steps 
to achieve the same accuracy compared to other methods. 

In this work we use method of lines coupled with RBFs to 
study the numerical solution of Coupled Drinfeld's-Sokolov- 
Wilson System: 

u, +3vv x = 0, 
v t + 2 v xxx + 2u v x + u x v = 0, (1.1) 

subject to following initial and boundary conditions 

u(x,0) = h(x), v(x,0) = g(x), a <x <b, (1.2) 

u(a, t) = h(t), u(b, t) = h 2 (t), v(a, t) = g\{t), v(b, t) = g 2 (t), 

(1-3) 

inspired by a recent approach of meshless MOL [15]. We gen- 
eralize this method to a system of coupled nonlinear PDEs. 

The method of lines (MOL) [? ] is generally recognized as 
a comprehensive and powerful approach to the numerical solu- 
tion of time-dependent PDEs. This method comprised of two 
steps: first, approximating the spatial derivatives; Second, then 
resulting system of semi- discrete ordinary differential equa- 
tions (ODEs) is integrated in time. Hence the method of lines 
approximates the solution of PDEs using ODEs integrators. 
In this paper, we will use the radial basis functions combined 
with the MOL to solve the system ( 11.11 ) inspired by [15]. As 
it evident from our results that this method possesses high ac- 
curacy and ease of implementation. The computed results are 
compared with the analytic solutions and good agreement is in- 
dicated. 

We find the numerical solution using radial basis functions 
(RBFs), particularly Hardy's multiquadric (MQ) and Guassian 
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(GA), as a spatial approximation for u, v and their derivatives 
in the the Coupled Drinfeld's-Sokolov-Wilson Equation. This 
transforms the system of partial differential equations to a sys- 
tem of first order differential equations in time. The solution can 
be obtained using a high order time integration scheme like the 
fourth order Runge Kutta method. Numerical results indicate 
that this method offers a highly accurate approximate solution, 
and it is easy to implement. This method does not require the 
grid generation as in the finite difference method, and the com- 
putation domain is composed of scattered collocation points. 
Higher order derivatives of u and v can be computed using the 
derivatives of infinitely continuously differentiable radial basis 
functions. 

This paper is organized as follows. In Section II we in- 
troduce the MOL combined with RBF method (MOL-RBF) 
for the Coupled Drinfeld's-Sokolov-Wilson Equation. Section 
III presents the implementation scheme and numerical experi- 
ments to demonstrate the effectiveness of the method proposed 
to solve the Coupled Drinfeld's-Sokolov-Wilson equation. Nu- 
merical results and comparison are given in Section IV and we 
conclude our work in Section V. 

2. Formulation of MOL-RBF method 

Radial basis function is a kind of function with the indepen- 
dent variable r, = r(x, Xi) = \\x - Some of the commonly 
used RBFs in the literature are: 

<pi(x) = (c 2 + r 2 ) 1/2 , Multiquadratics (MQ), 

<pi(x) = (c 2 + r 2 y 1/2 Inverse Multiquadratics (IMQ), 

(pi(x) - e~ cr > Gaussian (GA), 

where free parameter c is called the shape parameter of RBF. In 
the above definition x = (x, y) are the cartesian coordinates in 
£2 c R 2 , and the radius is given by 



and the matrix 



x - x. 



{(x-Xj) 2 + (y-yj) 2 }K 



where (xj, y j) is called the j source point of the RBF and is 
denoted by xj. We choose N nodes (x,-, i — 1 , 2, . . . , N) in P s c 
Q (J <9Q. Any given smooth function can be represented as a 
linear combination of RBFs: 



u»{x i ) = Y j A i <f> i = <i> T {x)A, 



(2.1) 



where 



<S>(x) = [(p 1 (x),<h.(x),...,<[> N (x)]' 1 
A - [A\, Ai, ■ ■ ■ , A N ] T . 



Which can be written as 



Aa = u, 



where 





<f> T ( Xl ) 




<Pi(xi) 


02(*l) 






A = 


4> r (x 2 ) 




(pi(x 2 ) 


<Pl(X2) 


■ ■ ■ 4>n(X2) 






. <I> T (x N ) . 




. <t>\(x N ) 


<p2(x N ) 


■ ■ ■ <Pn(xn) - 




is called the interpolation matrix, 


consisting 


of functions form 



ing the basis of the approximation space. It follows from 
Eq. d2.U and Aa = u that 



where 



u N (x) = <$> T (x) A~ x u = V(x) u, 



V(x) = <^ T (x)A- 1 . 



The convergence of RBF interpolation is given by the theo- 
rems in 1 19l 20ll : 



Assuming {xj} N , are source points in f s which is convex, 



the radial distance is defined as 



6 := <5(Q, P s ) = max min ||x - Xj\\2, 

xeCl l<i<N 



we have 



\\u N {x) 



u(x)\\ * Otf 15 ) 



(2.2) 



(2.3) 



U = [U\, tl2, 



,U N ] 



where < 77 < 1 is a real number and 77 = exp(-6>) with 9 > 0. 
From (12.3b it is clear that the parameter c and radial distance S 
affect the rate of convergence. 

The exponential convergence proofs in applying RBFs in 
Sobolov space was given by Yoon lfl6ll . spectral convergence 
of the method in the limit of flat RBFs was given by Fornberg 
et al.[2]. The exponential convergence rate was verified numer- 
ically by Fedseyev et al. iTTtH . The exponential convergence 
cited above is limited to certain classes of functions that are 
smooth enough and well-behaved in the domain of approxima- 
tion. 

In 1971, Hardy [5] developed multi-quadric MQ to approx- 
imate two-dimensional geographical surfaces. In Franke's ^ 
review paper, the MQ was rated one of the best methods among 
29 scattered data interpolation schemes based on their accuracy, 
stability, efficiency, ease of implementation, and memory re- 
quirement. Further, the interpolation matrix for MQ is invert- 
ible. In 1990, since Kansa [? ? ] modified the MQ for the 
solution of elliptic, parabolic and hyperbolic type PDEs, radial 
basis functions has been used to solve partial differential equa- 
tions numerically JHS0,[l3, 15]. The accuracy of MQ depends 
on the choice of a user defined parameter c called the shape pa- 
rameter that affects the shape of the RBFs. Golberg, Chen, and 
Karur [4] and Hickernell and Hon applied the technique of 
cross validation to obtain an optimal value of the shape param- 
eter c. 

The non-singularity of the collocation matrix A depends on 
the properties of RBFs used. According to IU8I1 . the matrix 
A is conditionally positive definite for MQ radial basis func- 
tions. This fact guarantees the non-singularity of the matrix A 
for distinct supporting points. This section includes details of 
the semi-discrete method, a method of line combined with the 
RBFs as in recent work fl5l, l7ll. 
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First, we approximate the space derivatives using radial ba- 
sis functions. The domain Q is divided into collocation points 
Xi, i = l,2,,..,N. Out of these points x,, i = 1,N belongs to 
the boundary <9Q. 
Let 



«"(**) =Z% 1 X}<Krij) = KxyA.u 
v N {xi) =Z% 1 A n j <f>(r iJ ) = 4>(x) T A 2 , 

where ry = \\xj - xj\\. In the matrix form 

u = A At, v = A Az, 
where u = [hi, u 2 , . ■ . , un], and v = [vi, v 2 , ■ • ■ , vjy] and 



(2.4) 
(2.5) 

(2.6) 





<f> T (xi) 






4> 2 (xi) . . 




A = 


<f> T (x 2 ) 




<pi(x 2 ) 


<f> 2 (x 2 ) ■ • 


• <Pn{x 2 ) 




. V(x N ) . 




. (pi(x N ) 


<p 2 (x N ) . . 


■ 4>n(xn) 



Therefore, from (12.4b and (12.6b . if the collocation matrix [A] is 
non-singular, it follows that 

u N (x) = (f> T (rj)A- l u = W(x)u, v N (x) = r (r/M -1 « = W(x)v, 

(2.7) 

where W(x) = <p T {r)A~ l . The non-singularity of the collocation 
matrix [A] depends on the properties of RBFs used. According 
to IHi . the matrix [A] is conditionally positive definite for MQ 
RBFs. This fact guarantees the non-singularity of the matrix 
[A] for distinct supporting points. 

Applying ( 12.7b on ( 11. lb . and collocation on each node xu we 
obtain 



dt 
dVj 
dt 



+ 3vi(W x (Xi)v) = 0, (2.8) 
+ 2(W xxx ( Xi )v) + 2u{W x ( Xi )v) + (W x {xdu)v = 0, (2.9) 
where M,(f) = u t , v,(f) = v,- and 

WAxd = [W u (xd Wuixt) . . . W lx (xi)l 
d 

WjAxi) = — U' ; U ; ). 

W xxx {Xi) = [Wux x (Xi) WuxxiXi) ■ ■ ■ Wuxx(Xi)], 

W jxxx (xd = g^WjiXi), 
for /' = 1 , 2, . . . , N. Then we rewrite the system ( 12.8b as 



dU 

— +3V* (W X V) = 0, 
dt 

dV 



(2.10) 



— + 2(W XXX V) + 2U * (W X V) + (W X U) * V = 0, (2.1 1) 
dt 

where * shows the component by component multiplicaton. 
Equivalently, this system can be written as 



dU dV 
— =F 1 (U), —=F 2 (V), 
dt dt 



where 
and 



FAU) = -3V * (W X V) 



F 2 (V) = - (2(W XXX V) + 2U * (W X V) + (W X U) * V) . 
The initial and boundary conditions translate to: 

U[t ] = [k°(jci), u°(x 2 ), u°(x N )] T (2.13) 



V(t ) = [v°( Xl ), v°(x 2 ), v»(x N )Y (2.14) 



0/„ si!" 



and 



(2.12) 



«i(0 = hit), u N (f) = h 2 (t), vi(0 = gi{t), v N (t) = g 2 (t). 

(2.15) 

To solve Eq. (12.12b - d2.15b . the fourth order Runge-Kutta 
scheme is applied as: 

K n =F l (U n ,t n ), K n = Fi{u n +*K u 

K 13 = F] (u n + jK 12 ) , K U = Fi (U" + AtK 13 ) 

K 2 i = F 2 (V, t„) , K 22 = F 2 Jy" + y K 2l J (2.16) 

K 23 = F 2 (v" + j ^22 J , K 24 = F 2 (V" + AtK 23 ) 

v n + l =U n + _± (Kn+ + 2Ki3 + Ru) 



yn+1 = yn + J_ {Kii + 2Ki2 + + ^ 

6 

Then, with the proper initial and boundary conditions, the com- 
putation can be carried out step by step. In the next section we 
implement our scheme to compute the solution of two examples 
of nonlinear coupled partial differential equations. 



3. Numerical Experiments 

In this section we apply the MOL-RBF method to the Cou- 
pled Drinfeld's-Sokolov-Wilson equations [14] and compare 
our results with other work in literature. 

Example . Consider the Coupled Drinfeld's-Sokolov-Wilson 
Equation with the following initial 

u(x,Q) = (c-4k)/2 + 3k 2 sech 2 (kx), v(x,0) = 2k Vc/2 sech(fex), 
and boundary conditions 

u{a, t) = (c - 4A:)/2 + 3 k 2 sech 2 (£(a - ct)), 

u(b, t) = (c - 4&)/2 + 3k 2 sech 2 (k(b - ct)), 

v(a, t) = 2k Vc/2 sech(Ai(a - ct)), 

v(b, t) = 2k 4cj2 sech(k(b - ct)) . 

The analytical solution |8j] is 

u(x,t) = (c-4k)/2 + 3k 2 sech 2 (k(x-ct)), 
v(x, t) = 2k s/c/2 sech(£(x - ct)) 



Numerical simulations are carried out on the interval [-4,4] 
from t = to t = 0.5 with step size h = 0.1 and time step 
0.0001, using Gaussian, Multiquadric, and Inverse Multi- 
quadric RBFs. Initial and boundary conditions are obtained 
from the exact solution in 111 411 . Table Q3-|4] shows the error for 
different values of c and k for GA and MQ. 



MoL-RBF Error u(x,t) 



-4-3-2-101234 



Figure 1 : Error with GA 



MoL-RBF Error v(x,t) 



-4-3-2-101234 



Figure 2: Error with GA 



4. Results and discussion 

We have applied meshless MoL to splve the coupled 
Drinfeld's-Sokolov-Wilson system using MoL-RBF and results 
are compared in TablesQJUfor different values of c and k. Ac- 
curacy was compared in terms of L M and RMS errors. As seen 
in Table[T]GA has accuracy higher than MQ. Since the accuracy 
depends upon the number of nodes and the value of the shape 
parameter c and we have obtained these values numerically by 
observing condition number of the matrix A and the error by 
keeping the number of nodes fixed. 



5. Conclusion 

In this work we constructed a MoL scheme combined with 
RBFs to solve system of two nonlinear partial differential equa- 
tions. This system has been solved in the literature using 
the Variational Iteration Method (VIM), Homotopy Analysis 
Method (HAM) and Homotopy Perturbation Methods (HPM). 
We have compared our scheme with the exact solutions taken 
from the cited literature. 
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Table 1 : Error comparison u(x, i) with the exact solution for c = 



1, k = 



0.001 





tlx 


0.1 


0.2 




0.3 




0.4 


0.5 




0.1 


4.5230 X 10~ 13 


8.6042 X 10" 


-15 


1.3291 X 10" 




3.5669 X 10" 12 


6.69920 X 10" 12 


MQ 


0.2 


1.3524 X 10" 12 


1.8085 x 10" 




1.3696 x 10" 




3.2862 X 10" ]4 


2.19896 x 10" 12 


C = 0.0001 


0.3 


2.2525 X 10 - ' 2 


3.6090 x 10" 


■12 


4.0690 X 10" 


■12 


3.6329 x lO - -' 2 


2.3009 X 10 - ' 2 




0.4 


3.1523 x 10~ 12 


5.4087 X 10" 


■12 


6.7682 X 10" 


■12 


7.2331 x 10~ 12 


6.8008 X 10~ 12 




0.5 


4.0521 x 10 - ' 2 


7.2083 X 10" 


■12 


9.4666 X 10" 


■12 


1.0832 x lO - ^ ' 


1.1300 X 10 - ' ' 




0.1 


8.10463 X 10~ 14 


2.2204 X 10" 




4.2293 X 10" 




6.84061 X 10~ 13 


1.0050 X 10~ 12 




0.2 


1.3206 X 10~ 13 


3.2413 x 10" 




5.7592 X 10" 




8.8817X 10~ 13 


1.2602 x 10~ 12 


IMQ 


0.3 


1.8307 X 10 -13 


4.2615 X 10" 


■13 


7.2891 X 10" 


■13 


1.0922 X 10 -12 


1.5153 X 10 -12 


C = 1 x 10 -7 


0.4 


2.8960 x 10" 13 


6.3921 X 10" 


■13 


1.0483 x 10" 


■12 


1.5184 x 10" 12 


2.0480 x 10" 12 




0.5 


3.4061 x 10 -13 


7.4124 X 10" 


■13 


1.2012 x 10" 


■12 


1.7225 X 10" i2 


2.3031 x 10" 12 




0.1 


2.9976 X 10 -14 







8.9983 x 10" 


-14 


2.3997 X 10 -13 


4.4997 X 10 -13 




0.2 


8.9983 X 10 -14 


1.1996X 10" 


■13 


8.9983 x 10" 


■14 





1.4999 X 10" 13 


GA 


0.3 


1.4999 X 10" 13 


2.3997 X 10" 


■13 


2.6995 x 10" 


■13 


2.3997 X 10" 13 


1.4999 X 10" 13 


C = 3400 


0.4 


2.0999 X 10" 13 


3.5999 X 10" 


13 


4.4997 X 10" 


■13 


4.7994 X 10" 13 


4.4997 X 10" 13 




0.5 


2.7000 x 10" 13 


4.8000 X 10" 


■13 


6.2999 x 10" 


■13 


7.1998 X 10 -13 


7.4995 X 10 -13 



Table 2: Error comparison for Drenfeld system with the exact solution for c = 4, k = 0.01 







0.1 


0.2 


0.3 


0.4 


0.5 




0.1 


3.5185 X 10 -8 


2.1287 x 10" 7 


5.3262 X 10" 7 


9.9565 X 10 -7 


1.6006 xlO -6 


MQ 


0.2 


8.1456 X 10 -10 


1.4087 x 10~ 7 


4.2467 X 10" 7 


8.5168 xlO -7 


1.4207 x 10" 6 


C = 0.0001 


0.3 


3.6814 X 10~ 8 


6.8879 x 10~ 8 


3.1672 X 10 -7 


7.0771 x 10" 7 


1.2407 x 10" 6 




0.4 


7.2813 X 10" 8 


3.1191 x I0" 9 


2.0877 Xl0 -7 


5.6374 X 10" 7 


1.0608 xlO -6 




0.5 


1.0881 X 10 -7 


7.5117 x 10 -8 


1.0081 X 10" 7 


4.1976 X 10" 7 


8.8087 X 10" 7 




0.1 


2.4470 X 10" 9 


1.4493 X 10" 8 


3.6138 XIO" 8 


6.7378 x 10" 8 


1.0820 X 10" 7 




0.2 


9.4369 X 10" 11 


9.7884 x 10 -9 


2.9081 X 10" 8 


5.7969 x 10" 8 


9.6450 X 10" 8 


IMQ 


0.3 


2.2578 x 10~ 9 


5.0841 x 10" 9 


2.2024 X 10 -8 


4.8562 XlO -8 


8.4693 x 10~ 8 


C = 10" 7 


0.4 


4.6096 x 10~ 9 


3.8058 x 10" 10 


1.4969 xlO -8 


3.9157 X 10 -8 


7.2938 x 10 -8 




0.5 


6.9607 x 10 -9 


4.3216 x 10" 9 


7.9168 x 10" 9 


2.9753 XlO" 8 


6.1185 xlO -8 




0.1 


2.3999 x 10" 9 


1.4399 x 10" 8 


3.5997 xlO" 8 


6.7189 X 10" 8 


1.0797 X 10 -7 




0.2 





9.5997 X 10" 9 


2.8798 x 10 -8 


5.7592 X 10" 8 


9.5979 x 10' 8 


GA 


0.3 


2.3999 x 10 -9 


4.7999 X !0" 9 


2.1598 xlO" 8 


4.7994 x 10 -8 


8.3983 x 10" 8 


C = 2700 


0.4 


4.7999 x 10~ 9 





1.4399 x 10"" 8 


3.8395 XlO -8 


7.1986X 10 -8 




0.5 


7.1998 x 10 -9 


4.7999 X I0 -9 


7.1996 xlO -9 


2.8797 X 10" 8 


5.999 x 10 -8 



Although the optimal value of shape parameter is open prob- 
lem, this proposed method has many benefits. First, it has 
higher accuracy for GA RBFs. Second, it uses the meshless 
property of RBFs, as the mesh generation is expensive for many 
complicated geometries and for when there are singularities. 
Third, for system of coupled PDEs, this proposed method has 
accuracy comparable to other methods yet it is easier to imple- 
ment for higher order PDEs due to infinitely differentiability of 
RBFs used. 

Further investigations can be carried out for larger number 
of nodes, optimal value of shape parameter, using adaptive step 
Runge-Kutta method or other multistep time integrators. Re- 
cently some work has been done on adaptive RBFs methods. 
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Table 3: v(x, i) Error comparison for Drenfeld system with the exact solution for c = 0.1, k = 0.001 





tlx 


0.1 




0.2 




0.3 




0.4 




0.5 






0.1 


6.5428 X 10" 


14 


1.1693 X 10" 


13 


1.5573 X 10" 


13 


1.7871 X 10" 


13 


1.8790 X 10" 


13 


MQ 


0.2 


1.3783 X 10" 


13 


2.6174 X 10" 


13 


3.7467 x 10" 


13 


4.6779 X 10" 


13 


5.4992 x 10" 


13 


C = 0.0001 


0.3 


2.1028 x 10" 


13 


4.0664 x 10" 


13 


5.9372 X 10" 


13 


7.5756 X 10" 


13 


9.1212 x 10" 


13 




0.4 


2.8279 X 10" 


13 


5.5165 x 10" 


13 


8.1294 X 10" 


13 


1.0475 X 10" 


12 


1.2746 x 10" 


12 




0.5 


3.5539 X 10" 


13 


6.9683 x 10" 


13 


1.0324 X 10" 


13 


1.3378 X 10" 


12 


1.6374 X 10" 


12 




0.1 


2.5528 X 10" 


13 


6.8342 x 10" 


13 


1.2824 X 10" 


12 


2.0582 x 10" 


13 


3.0049 x 10" 


12 




0.2 


4.4478 x 10" 


13 


1.0623 x 10" 


12 


1.8487 X 10" 


12 


2.8156 x 10" 


12 


3.9515 x 10" 


12 


IMQ 


0.3 


6.7614 x 10" 


13 


1.5248 X 10" 


12 


2.5402 x 10" 


12 


3.7401 X 10" 


12 


5.1066 x 10" 


12 


C = 10 -7 


0.4 


9.7163 x 10" 


13 


2.1 156 X 10" 


12 


3.4239 x 10" 


12 


4.9207 X 10" 


12 


6.5817 X 10" 


12 




0.5 


1.3551 x 10" 


12 


2.8823 x 10" 


12 


4.5714 x 10" 


12 


6.4529 x 10" 


12 


8.4963 x 10" 


12 




0.1 


1.6994 x 10" 


12 


3.2199 x 10" 


12 


4.5615x10" 


12 


5.7243 x 10" 


12 


6.7082 x 10" 


12 




0.2 


3.4882 x 10" 


12 


6.7976 X 10" 


12 


9.9281 x 10" 


12 


1.2879 X 10" 


11 


1.5652 x 10" 


11 


GA 


0.3 


5.2771 x 10" 


12 


1.0375 X 10" 


11 


1.5294 x 10" 


11 


2.0035 X 10" 


11 


2.4596 X 10" 


11 


C = 2700 


0.4 


7.0659 x 10" 


12 


1.3953 X 10" 


11 


2.0661 x 10" 


11 


2.7190 X 10" 


11 


3.3541 X 10" 


11 




0.5 


8.8548 x 10" 


12 


5.8014 x 10" 


11 


2.6027 X 10" 


11 


3.4346 X 10" 


11 


4.2485 X 10" 


11 



MQ 

C = 0.0001 



OA 

C = 2700 



Table 4: v(x, t) Error comparison for Drenfeld system with the exact solution for c = 0.0001, k = 0.01 



o.i 



0.2 



0.3 



0.4 



0.4 
0.5 



2.2379 X 10~ n 
4.4328 X 10~ n 
6.6276 X 10~ 11 
8.8225 X 10~ n 
1.1017 X 10" 10 



4.5619 X 10" 
8.9517 X 10" 
1.3341 X 10" 
1.7731 X 10" 
2.2120 X 10" 



6.9696 X 10" 
1.3552 X 10" 
2.0134 X 10" 
2.67 16 X 10" 
3.3298 X 10" 



9.4682 X 10" 
1.8247 X 10" 
2.7027 X 10" 
3.5806 X 10" 
4.4585 X 10" 



2.0340 X 10 



1.0306 X 10 



"12 



II 



4.0619 X 10" 
8.1481 X 10" 
1.2258 X 10" 
1.6406X 10- 
2.0605 X 10" 



1 1 



6.0816 X 10" 
1.2208 X 10" 
1.8372 X 10" 
2.4591 X 10" 
3.0888 X 10" 



1 1 



8.0993 X 10- 
1.6271 X 10- 
2.4492 X 10 - 
3.2787 X 10- 
4.1185x10- 



14 



0.3 
0.4 
0.5 



5.6457 X 10" 
1.1291 X 10- 13 
1.6937 X 10" 13 



13 



1.1290X 10" 
2.2582X 10~ 13 
3.3874X 10 13 



13 



1.6935 X 10 
3.3873 X 10- 13 
5.0811 X 10- 13 



2.2583 X 10 



■13 



4.5166 X 10 



13 



6.7748 X 10 



13 



2.2579 X 10- 
4.5163 X 10- 
6.7747 X 10- 
9.0330 X 10- 
1.1291 X 10" 



0.5 



1.2050 X 10- 
2.3025 X 10- 
3.3999 X 10- 
4.4973 X 10- 
5.5947 X 10" 



1.0108 X 10- 
2.0324 X 10- 
3.0600 X 10" 
4.0969 X 10- 
5.1466 X 10- 



2.8223 X 10- 
5.6452 X 10 - 
8.4682 X 10- 
1.1291 X 10- 
1.4114X 10- 
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Figure 4: c = 0.0001 



Figure 6: c = 4 
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Figure 7: u(x, t) with GA 
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Figure 8: v(x, t) with GA 




Figure 9: u(x, t) with MQ 

X 10"" 




Figure 10: v(x, i) with MQ 
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Figure 11: u(x, t) with IMQ 

x10~ 4 
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Figure 12: v(x, f) with IMQ 
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Figure 13: Error vs. Shape Parameter 
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Figure 14: Error in u with GA 
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Figure 15: Optimal Shape for GA, k = 
0.01, c = 0.0001 




Figure 16: Optimal Shape for MQ, k = 
0.01, c = 0.0001 
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